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Abstract 



A new mathematical model for the dynamics of prion proliferation involving an ordinary differential equation 
coupled with a partial integro-differential equation is analyzed, continuing the work in We show the well- 
posedness of this problem in its natural phase space := X ((xq, oo); xdx), i.e. there is a unique global 
semiflow on Z+ associated to the problem. 

A theorem of threshold type is derived for this model which is typical for mathematical epidemics. If a certain 
combination of kinetic parameters is below or at the threshold, there is a unique steady state, the disease-free 
equilibrium, which is globally asymptotically stable in Z+; above the threshold it is unstable, and there is another 
unique steady state, the disease equilibrium, which inherits that property. 



This paper was initiated while the second author was visiting the Department of Mathematics, Vanderbilt 
University, Nashville, Tennessee in 2003/04. He wants to express his thanks to the department for kind hospitality 
and for financial support. 



In this paper we continue our analysis, begun in [9], of a recent model describing the proliferation 
of prions. This model has been introduced in Greer, Pujo-Menjouet and Webb jjj, based on the 
works of Masel, Jansen and Nowak [SJ, Nowak, Krakauer, Klug and May |Zj and others. For 
comprehensive explanations and discussions of the model and the relevant biochemical literature 
we refer to 0]. Here we only give a very short description of the model. 

Prions are proteins that are believed to be responsible for certain diseases like BSE and the 
Creutzfeld-Jacob disease. There are two basic forms of prions of interest here, the Prion Protein 
Cellular PrP c and the Prion Protein Scrapie PrP Sc . The single molecule proteins PrP c , also 
called monomers in the sequel, are protease resistent proteins which have a cell protective function 
and are produced by the body, regularly. On the other hand, the infectious prion PrP Sc is a 
string- like polymer formed of monomeric PrP c . Above a critical chain length xq > the polymers 
are more stable than the PrP c , and they can grow to chains containing thousands of monomers. 
p r pSc ^he ability to replicate by splitting, we assume binary splitting here. 

So there are three main processes which govern the dynamics of prions in this model. 

• growth in length by polymerization with rate r > 0; 

• binary splitting with rate (3{x) > 0, a polymer of length x > splits into one of length 
< y < x and one of length x — y with probability n(y, x)\ 

• natural degradation with rate 7 > for the monomers and with rate fi(x) for the polymers 
with length x. 

The model proposed in |Jj further assumes that polymers of length < x < xq immediately 
decompose completely into monomers. This reflects the assumption that PrP Sc polymers are 
unbranched and form a simple a-helix with xq monomer units per turn. An a-helix of length less 
than xq is incomplete and thus is much less stable. Denoting the numbers of monomers at time t 
by V(t) and the density of polymers by u(t,x), we obtain the following model equations. 
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1. Introduction and Main Results 




u(t, x)dx + 2 




(1.1) 



d t u{t, x) + TV(t)d x u(t 7 x) + (n(x) + P(x))u(t, x) = 2 




V(0) = V" o >0, u(t,x ) = 0, u{0,x) = u {x), 
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where t > and xq < x < oo. Here A > is a constant background source of monomers. Observe 
that the splitting function n{y, x) should satisfy the following properties. 



K(y,x)>0, K(y,x) = k(x - y,x), I n{y,x)dy = l, 

Jq 

for all x > xq, y > 0, and n(y, x) = if y > x or x < xq. Note that these conditions imply 

2 / yn(y,x)dy — x, x > 0. 
Jo 

In fact, 

f-X f-X f-X 

yn(y,x)dy= / yn(y,x)dy+ / yn(x - y,x)dy 
Jo Jo 

X f-X r-X 

yn(y,x)dy+ (x - y)n(y 1 x)dy = a; / n(y,x)dy = x. 



This implies that mass does not change via the splitting process, and by a simple computation 
we obtain the following relation for the total number of monomers in the system. 

-r[^)+/ xu(t,x)dx] = A - 7V(f) - / x(i(x)u(t,x)dx, t>0. 

®t J x J Xq 

In it is further assumed that splitting is equi-distributed (polymer chains are equally likely to 
split at all locations), and that the rate of splitting is proportional to length. This reflects again the 
hypothesis that polymers form a-helices and are not folded in more complicated configurations, 
which would make certain segments of the chain less likely to split than others. Therefore, we 
make the further assumptions 

ft(y, x) = 1/x if x > xq and < y < x, k(jj, x) — elsewhere , 

P{x) = fix is linear, and fi{x) = [i constant. Then the model contains only 6 parameters, and 
can even be reduced to a system of 3 ordinary differential equations. In fact, introduce the new 
functions 



U(t)= / u(t,y)dy and P(t) = / yu{t,y)dy, 

Jxq J Xq 

representing the total number of polymers, and the total number of monomers in polymers at time 
t, respectively. Integrating the equation for u(t,x) over [xo,oo) we get 

d f°° f°° 

-U{t) = -rV{t)u{t,x)\Z-lJLU(t)-pP{t) + W / u(t,y)dydx 



dt 

/>OC 

= -tiU(t)-0P(t)+20 u(t,y)(y-x )dy 

J x 

= -iiU(t) - /3P(t) + 2(3P(t) - 2/3x U(t), 



hence 



U(t) = -(ji + 2/3x )U(t) + f3P{t). 
Multiplying the equation for u(t, x) by x, integration yields 
d 



di P(t) = -TV(t)(xu(t,x)\™- I u(t,y)dy) 



-fiP(t)-P u{t, x)x z dx + 2/3 / x u{t,y)dydx 

•J XQ J Xq J X 

= TV(t)U(t) - fiP(t) - pi u(t,x)x 2 dx + (3 u{t,y)(y 2 - x 2 )dy 



hence 



rV(t)U(t) - (iP(t) - (3xlU{t), 

P{t) = rU(t)V(t) - fiP(t) - pxlU{t). 
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Thus we obtain the following closed model involving only ordinary differential equations. 

U = 0P-I1U- 2f3x U 
(1.2) V = A - ~fV - tUV + (3x 2 U 

P = rUV-fiP-PxlU 

with initial conditions 

U(0) = U Q > 0, V(0) = V > 0, P(0) - P > x U . 

This way the partial differential equation for the density u(t, x) decouples from the ordinary dif- 
ferential equations. Once the solutions of p. 2(1 are known, one has to solve only a linear partial 
integro-differential eqution to obtain u(t, x). The system (1.2) is identical to the "basic virus 
dynamics model" that is discussed at length in [Hj. 

Concerning the ode-system ((1.2(1 we have the following result from Priiss, Pujo-Menjouet, Webb 
and Zacher 

Theorem 1.1. Suppose xo, /3,7, A, /x, r > are given constants. Then the system induces a 

global semiflow on the set K = {(U, V, P) £ I 3 : U, V, P — xqU > 0}. There is precisely one disease 
free equilibrium (0, A/7,0) which is globally exponentially stable if and only if /i + Xo/3 > y A/Jr/7, 
and asymptotically stable in case of equality. On the other hand, if fi + xq/3 < \J A/3t/7 there is 
the unique disease equilibrium 

/ \pT-j(p + f3x ) 2 {n + f3x )) 2 A/3T- 7 (/i + /3a:o) 2 \ 
V ht(h + 2(3x ) ' (3t ' far ) 

which is globally exponentially stable in K \ [{0} x R + x {0}] . 

ft is the purpose of this paper to study the full system p. 1(1 under the assumptions of equi- 
distributed splitting, linear splitting rate, and constant rates of degradation. 

Since V(t) + xu(t, x)dx is the total number of monomers in the system, which should be finite 
at any time, it seems reasonable to study 1(1.1(1 in the standard cone Z + := R + x L^((xo, 00); xdx) 
of the Banach space Z :=Rx Li((xo,oo);xdx). The following theorem summarizes our results. 

Theorem 1.2. Assume equi- distributed splitting with linear splitting rate P(x) = fix and constant 
degradation rates 7 and ji{x) = /i. Suppose A, r, /3, 7, /x, Xq > 0. Then ((1.1(1 generates a global 
semiflow in the natural phase space Z + . Furthermore, 

(i) if A/3t/ 7 < (/i + Pxq) 2 , then the disease-free equilibrium z — (A/7,0) is globally asymptotically 
stable in Z + , and even exponentially in the case of strict inequality; 

(ii) if A^r/j > (/i + f3xo) 2 , then there is a unique disease equilibrium — (T4,u») which is 
globally asymptotically stable in Z + \ (R + x {0}). It is given by 

v _ (m + /?^o) 2 u _ W A/3r - 7 ( M + f3x Q ) 2 ^, P{x- x ) , 
(3t ' * [it [p, + f3xo){n + 2/3xo ) (j, + [3xq 
where $(r) = (r + r 2 /2) exp(— (r + r 2 /2)). 

The remaining part of this paper deals with the proof of this result. Recall that the function 
u(t) := rV(t) can be considered as known, by Theorem 1.1, and uj(t) — > ujoo exponentially, where 
either = A/7 in the disease-free or ojoc = (/Lt + (3xo) 2 / (3 in the disease case. Hence we have to 
solve a linear nonautonomous partial integro-differential equation of first order. For this we shall 
use standard techniques from the theory of Co-semigroups and we refer to the monograph Arendt, 
Batty, Hieber and Neubrander [2] as a general reference for the results employed below. 

We proceed in four steps. First we study the autonomous case where to = lo^. In Section 2 
we show that there is a unique Co-semigroup T(t) — e~ Lt associated with the pde-part of 1(1.1(1 
in X — L\((xq, 00); xdx), which is positive and contractive, and even exponentially stable in the 
disease-free case. The resolvent of L is shown to be compact in Section 3, hence L has only point 
spectrum in the closed right half-plane. In the disease case, we further show that is the only 
eigenvalue of L on the imaginary axis, it is simple and so the ergodic projection V onto the kernel 
N(L) of L along the range R{L) of L exists and is rank one. We compute an element e G N(L) 
which is positive. A result of Arendt, Batty, Lubich and Phong |2] then shows that T(t) is strongly 
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ergodic, i.e. Hindoo T(t) — V strongly in X. Wellposedness of the nonautonomous problem is 
proved in Section 4 by means of monotone convergence, it is shown that the evolution operator 
exists and is bounded. Moreover, bounds for d x u(t, •) in X are derived. Finally, in Section 5 we 
put together these results to prove Theorem 1.2. 

While we assume throughout that f3(x) = fix, /i(x) = [i (constant), and yn{x,y) — 1 for x < 
y, y > xq, k(x, y) = elsewhere, our methods extend to versions of (1.1) where these assumptions 
do not hold. We do not carry out these generalizations since it is not clear which would be 
biologically reasonable. On the other hand, the equation discussed in this paper 

d t u(t, x) = ~TV(t)d x u(t, x) — (/i + (3x)u(t, x) + 2/3 / u(t, y)dy 

J X 

for x > xo, t > 0, with initial and boundary data as in (1.1), can be solved with an integral 
transformation followed by the method of characteristics. Namely, define 



v(t,x)= / / u(t,0dtdy= / (Z-x)u{t,Z)dt, d£v(t,x) = u{t,x). 

J x J y J x 

Then a computation shows that v solves the first order partial differential equation without integral 
term 

d t v(t, x) — —TV(t)d x v(t, x) — (fi + (3x)v(t, x) 

for x > xo, t > 0, with initial data v(Q,x) obtained by integrating uq twice and boundary data 
v(t,x ) = P(t) — xolf(t). The equation for v may be solved by the method of characteristics, and 
u is recovered from d x v{t,x) — u{t,x). The solution depends on the initial data in the region 
{(a;, t) | x > xq + r J Q V(s)ds } and on the boundary data in the complement of this region. Since 
V(t) always has a positive limit, it is evident that the contribution from the initial data is swept 
out towards large x-values and decays exponentially, in fact, at a rate like e for some e > 0. If 
the disease-free state is stable, then (P(t), U(t)) — > (0, 0) as t — > oo, which implies that the solution 
u converges to zero also in the region where it depends on the boundary data. In the case of a 
positive disease equilibrium, P(t) — xoU(t) has a positive limit as t — * oo, which determine the 
limiting equilibrium distribution given in Theorem 1.2. This method breaks down if /?(•), //(•), 
or k(-, ■) have more complicated forms, as the reader will readily confirm. 

2. The Linear Automomous Problem 
2.1. Functional Analytic Setting. We consider the problem 

poo 

(2.1) d t u(t, x) + Ldd x u(t, x) + (p + I3x)u(t, x) = 2/3 / u(t,y)dy, 

J X 

u(0,x) = uq(x), u(t,xa) — 0, t > 0, x > xo- 
Set w(t, x) = u(t, x + xq), x > 0. Then this problem becomes the following one on IR_|_. 

/>oo 

(2.2) d t w(t,x) +ujd x w(t,x) + (fi + f3x)w(t,x) =2(3 w(t,y)dy, 

J X 

w(0, x) = g{x) := uq{x + xq), w{t, 0) = 0, t > 0, x > 0. 
Here we have set /io = [i + /3xq. uj plays the role of tV at oo, i.e. 

lu = tV(oo) = At/7 

in the disease-free case or 

uj = tV{oo) = (/x + Px ) 2 /P = 

in the disease case. 

We want to study l|2.2|l in the basic space X = Li(R + ; (a + x)dx), where we choose as the norm 

||w|| = a|w|i + \xw\i, 

with a > to be determined later. We define two linear operators in X by means of 

Au(x) = wu'{x) + (/io + (3x)u{x), x € R+, 
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with domain 

D(A) = {ue Wt{R+)nX : x 2 u e L 1 (M + ),a;ti / (a;) € u(0) = 0}, 

and 

/•OO 

Bu(i) = 2/3 / u(y)dy, £>(5) = D(A). 



Both operators are well-defined and linear, B will be considered as a perturbation of A. 
2.2. m-Accretivity of A. We have 

I />oo 

sgn mix = u> \u\' dx + fio\u\i + (3\xu\i 
Jo 

= fio\u\i + 0\xu\t, 

and 

/>OC />oc 

/ Ausgnuxdx = uj \u\'xdx + hq\xu\i + f3\x 2 u\i 
Jo Jo 

= —u)\u\i + ^o|*w|i + P\x 2 u\i. 

Employing the bracket in L\ this implies 

[Au, u]+ > (a/io - oj)\u\i + (a/3 + /zo)|arat|i > r]\\u\\, 

for some 77 > provided (1q > Lu/a. Hence for such a, A is strictly accretive, in particular closable. 

Next we compute the resolvent of A. The equation (A + A)u = / is equivalent to solving the 
ode 

(2.3) Xu(x) + wu'(x) + {no + (3x)u(x) = f(x), x > 0, 

with initial condition u(0) = 0. Therefore we obtain 

u=(X + A)- 1 f(x) = - f exp-{(\ + no)(x-y)/LO + P(x 2 - y 2 )/2oj}f(y)dy. 



If / € Xi(M+) then on easily obtains the estimate 

< I/Ix/CA + mo). 

If also xf e L 1 (R + ) then 

\a?u(x)\ < - [\-^+^ x -y)^(x 2 -y 2 )e'^ x2 'y 2 ^ 2 ^\f{y)\dy 
w Jo 

+ - [ X ] ye- f3(x - y)2y/2 "y\f(y)\dy, 
u Jo 

hence 

, 9 , 1 w 2w , „ 1 w. „ 

A < -T-— 7T / 1 + - m x f i- 

w A + ^ pe wp 2 

This shows that x 2 it S hence m £ Li(M. + ), and then by equation (|2.3[1 also it' G Li(R+) 

as well as iu' G Li(R+), i.e. it s D(A). This shows that ^4 is m-accretive. 

As a consequence we note that —A generates a Co-semigroup in X which is also positive and 
strictly contractive, hence exponentially stable. 

2.3. Accretivity of A — B. We have 

/>oo />oo -t 

I / u(x)dx\i < \xu\i, \x / u(x)dx\i < — \x 2 u\i, 

o X «/ X 

and therefore 

/>oo 

/ (Au — Bu) sgn(u)dx > /io|u|i + /3|a:u|i — 2/3|xu|i, 
Jo 

as well as 

poo 

(Au — Bu) sgn(u)xdx > —ui\u\i + fio\xu\\. 
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This yields 

[(A - B)u, u}+ > (fi a - lj)\u\i + (mo - /3a)\xu\i > 0, 

for all u e D{A), provided fio a > w and yu > /3a. Such a choice of a > is possible if and only if 
the condition w/^o < Mo//3 is met, i.e. if and only if 

w < Mo//? 

holds true. Now in the disease-free case we have u = At/7, while in the disease case u) = /Uq//3; 
then a = Ho/(3. Thus A — B will be strictly accretive in the disease-free case while it will be 
accretive only in the disease case. In the first case, the decay rate can easily be estimated not to 
be smaller than /j,q — y/X^r/j. 

2.4. Density of the Range of A — B. Let / e Li(R + ; (a + x)dx) be given and assume / > 0. 
Set Mi = (1 + A)^ 1 f and define the sequence u n inductively by means of 

u n+ i = Ml + (1 + A^BUn. 

Then u\ > 0, and u 2 — U\ = (1 + ^4) _1 _Bit 1 > 0, hence by induction u n+ \ > u n pointwise, since B 
is positive. This shows that the sequence of functions u n is nonnegative and increasing pointwise. 
Moreover, 

uju' n + (1 + fi Q + (3x)u n = f + 2(3 / u n -!(y)dy < f + 2(3 / u n (y)dy, 

J x J X 

which implies 

(1 + /j )K|i + (3\xu n \! < l/li + 2(3\xu n \!, 

and 

-oj\u n \i + (1 + no^xunlx + /3 2 |a; 2 u„|i < \xf\i + /3|a; 2 u„|i. 

Choosing a as above this yields an a priori bound for the sequence (u n ) 

\\u n \\ = a\u n \i + \xu n \ < C\\f\\, 

and therefore we may conclude by the monotone convergence theorem u n — > as n — > oo. If 
in addition x 2 f E Li(R + ) then we obtain in a similar way boundedncss of x 2 u n in X. This 
implies (1 + A — B)u n = f + B(u n -i — u n ) — > / in X as n — > oo, hence G _D(A — _B) and 
Moo = (1 + A — B)^ 1 f . Since L\ = L\ — L\ we may conclude R(l + A — B) = X, i.e. the closure 
of A — B is m-accretive. 



Remark 2.1. The above proof shows that the resolvent of A — B is positive, hence the semigroup 
generated by this operator will be as as well. 

2.5. Irreducibility. Suppose / £ X is nonnegative and u solves 

wu' + (A + + @x)u = 1 + 2(3 I u(y)dy, x > 0, 

J X 

with initial value u(0) = 0. If / ^ then let xi := inf supp /. We have 

u{x) = - f cxp-[(\ + n Q )(x-y)/co + p(x 2 -y 2 )/2cj}[f(y) + Bu(y)}dy. 
^ Jo 

Since we already know u(x) > 0, this formula implies u(x) > for all x > X\. But then 
u{y)dy > for all x > 0, and so so u(x) > for all x > 0. This proves the irreducibility 
of the semigroup generated by A — B. 
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2.6. A — B is not Closed. Unfortunately, the sum A — B is not closed. We show this by the 
following example. 

Example 2.2. Set u — x/x 3 where x denotes a cut-off function which is on [0,1] and 1 on 
[2, oo). Then u,u'u,xu € Li(R+), but x 2 u ^ Li(R+), and u(0) = 0. On the other hand, 



f(x) := luu'(x) + (X + (i + I3x)u(x) -2/3 u{y)dy 

J X 

= ujx'/x 3 - Zlox/x* + (A + (io)x/x 3 + Px/x 2 - 2/3 / x (y)dy/y 3 



Since 



x(x)/x 2 - 2 J x (y)dy/y 3 = x(x)/x 2 + x (y)/y 2 \T - / x\y)dy/y 2 
x'(y)dy/y 2 , 



we obtain 

/>oo 

/ = lux'(x)/x 3 - 3lu X (x)/x 4 + (A + no)x(x)/x 3 - / x'(y)dy/y 2 . 
Obviously, / as well as xf belong to Li(R+), so A — B with domain D(A) is not closed. 
2.7. Summary. Let us summarize what we have shown so far. 

Theorem 2.3. Suppose 0u> < (1q. Then problem ]2.ty) is well-posed in X — Li(R+; (a + x)dx) and 
admits an associated Co-semigroup T(t) — e~ Lt which is positive. If a is chosen from the interval 
a e [u>/(io , no /(3] then T(t) is nonexpansive. 

In the strictly disease free case u> — Ar/7 < (1q/(3, the semigroup T(t) is exponentially stable 
with type u>o(T) < — (i + y X/3t/j < 0. 

3. Asymptotic Behavior of the Autonomous Problem 



3.1. Compactness. Set L — A — B. Since L is m-accretive in X = Li(M+;(a + x)dx), the 
spectrum a(L) is contained in the closed right halfplane. We want to show that the resolvent of L 
is compact. For this purpose we derive another representation of (A + L)^ 1 for A > 0. Let f E X 
and set u = (A + L)^ 1 f . Then we obtain 



and 



u = (A + A)-\f + (A + A)~ x Bu, 



{\ + A)- 1 Bu = 2(3{X + A)- 1 [ u{y)dy] 



= 2 JL f X e -(^oKx~y)/^ e -f3(x*~y*)/2u J[ I u{r)dr]dy 



20 

id 



u(r)[ / e -(A+Po)(«-v)/" e -/3(x 
Jo 

+ 2 JL f u {r)[f e-t x+ ^-vV"e-^ 2 -v 2 )/^dy]dr 
u Jo Jo 

poo 

= k x (x) / u(r)dr + + Ay 1 [k x u], 



where 



Note that 



kx{x ) = 2 A I e -^+*>)(*-v)/" e -P(* 2 -y 2 )/^ dy . 



< k x (x) < ^ f e-^ + ^ x -^ u dy < W 



A + (io ' 
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i.e. k\ £ L 00 (R + ). We thus have the identity 

/>oo 

u(x) — k\(x) / u(y)dy^(\ + A)- 1 f(x)+u(\ + A)- 1 [k x u}^:g(x), 

J X 

and u(0) = 0. We may solve this equation for u to the result 

u(x) = g(x) - k x {x) J cxp ^- J k\(r)drj g{y)dy + k\(x) exp ^- J k\(s)ds^ < q x \f >, 
where 

i POO />00 

This way we have the representation 

(3.1) (X + L)-\f = (l- J R A )(A + A)- 1 [l+^A(A + L)- 1 ]/ + fc A (x)cxp (- £ k x (s)ds^ < q x \f >, 



with 



(R\g)(x) = k\(x) J exp^- J k\(r)drj g(y)dy. 



Next D(A) embeds compactly into X, hence (A + A)^ 1 is compact. From boundedness of k\ we 
may then conclude that (A + L)^ 1 is compact, as soon as we know that the Volterra operator R\ 
is bounded in X. 

To prove the latter we estimate as follows 

pOO f-X / f-X \ 

||i? A .g|| = J (a + x)kx(x)\j exp ^- y k x {r)drjg{y)dy\dx 

pOO pOO / pX \ 

< y \g(y)\[J (a + x)kx(x) exp ^- J k x {r)drj dx]dy 

pOO pOO / pX \ 

= y \g(y)\[( a + y) + J exp \~J k *( r ) dr ) dx \ d y 

pOO 

< C x / |fl(tf)|(o + y)d|, = CA||fl||, 



as we show now. 



jfe^) = M f X e -(X+^)(x-y)/^ e ^(x 2 -V 2 )/2^ dy 

^ Jo 

> M / e -(*+»o)v/u e -l3xy/u d 



- 1 ^ _ e -(A+ M +/3x)x/u^ 



> 



A + /i + /3a; 

2/3 (A + + /9a;)a;/w 



A + /i + /3a; 1 + (A + fi a + (3x)x/u>) 
2(3x 



uj + (A + /i + /3a;)a; ' 
by the elementary inequality 1 — e~ x > a;/(l + a;). This implies 

f-X f-X 

I k x {r)dr > 2/3 / rdr/(w + (A + ^ + /3r)r 

f x 2/3r + A + Mo /"* dr 

= / ; 7T^ dr - + Mo) 

L w + (A + na)r + (ir 2 



uj + (A + /i )r + /3r 2 



. . uj + (A + /i ).x + /3a; 2 
w + (A + ^ )y + /3y 2 
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since the second integral is bounded. This estimate finally yields 

f 00 exp (- f kx (r)dr) d X < e- f° "+ ( ^f ^ a + »>" 

Jy V Jy J Jy + fl )x + (3 X 2 

This completes the proof of compactness of the resolvent of L. 

3.2. Ergodicity. Since the resolvent of L is compact we know that the spectrum of L consists 
only of eigenvalues of finite multiplicity, these are poles of the resolvent of L. By accretivity of 
L we have the inequality |(A + £) _1 |b(x) < 1/ReA, ReA > 0, hence the resolvent can only have 
poles of first order on the imaginary axis. This shows that all eigenvalues on the imaginary axis 
are semisimple. Compactness of the resolvent implies also that the range of A + L is closed, for 
each A e C. In particular, we have the direct sum decomposition X — N{L) © R(L), i.e. ergodicity 
in the sense of Abel. 

Now we concentrate on the disease equilibrium which means a = Mo//3 an d u> — fJ-o/0- A 
function e(x) belongs to the kernel of L if 

/>oo 

ue'(x) + (// + /3x)e(x) - 2(3 j e(y)dy = 0, x > 0, e(0) = 0, 

J X 

or equivalently 

e"(x) + — (1 + —x)e'{x) + 3^e(x) =0, x > 0, e(0) = 0. 
Mo Mo Mo 

The scaling e(x) — v(J3x/ reduces this problem to 

v"(z) + (1 + z)v'(z) + 3v(z) = 0, z > 0, v(0) = 0. 

By the initial condition v(0) = 0, this shows that the kernel of L can be only one-dimensional, and 
a simple computation yields that 

v{z) = {z + z 2 /2)e- {z+z2 ' 2 \ z>0, 

is a solution. Therefore N{L) = span{e}, with e(x) = {(3/ /j,o) 2 v((3x/ Mo) 5 an d another simple 
computation yields 

/•OO 

/ (a + x)e(x)dx = 1. 
^o 

Since L is Fredholm with index zero, the kernel N(L*) of the dual of L has also a one-dimensional 
kernel which are the constant functions. The ergodic projection V onto the kernel of L along the 
range of L is then given by 



(3.2) Vu(x) = [ (a + x)u{x)dx]e{x) =< u\e* > e(x), x > 0. 

^o 

Suppose there are no other eigenvalues of L on the imaginary axis. Then L* also has no other 
eigenvalues on the imaginary axis, and then by the theorem of Arendt, Batty, Lubich and Phong 
we may conclude that 

e~ Lt u — > Vu as t — > oo, for each u G X, 

i.e. the semigroup generated by — L is strongly ergodic. 

We show now that there are in fact no eigenvalues other than on the imaginary axis. Suppose 
on the contrary that 

/>OC 

ipu(x) + lou'(x) + (/io + (3x)u{x) — 2(3 j u(y)dy, x > 0, u(0) = 0, 

J X 

w^O. Multiplying this equation with u/\u\, taking real parts, and integrating over M + we obtain 

pOO f-X 

(3.3) MoMi +fi\xu\x = 2/3Re / u(x) u{y)/\u(y)\dydx <2(3\xu\ u 

Jo Jo 

and similarly, multiplying with xu(x)/\u(x)\ we get 

/>0O f-X 

(3.4) — lo\u\i + fia\xu\i + (3\x 2 u\i = 2/3Re / u(x) / yu(y)/\u(y)\dydx < (3\x 2 u\i. 

Jo Jo 
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Multiplying the first inequality with a = po/(3 and adding the second we arrive at a contradiction 
if at least one of the inequalities i|3.3|) . i|3.4f> is strict. Hence we must have 



fOO f>X 

Re / u(x) / u(y)/\u(y)\dydx=\xu\i, 



io Jo 
which implies with argu(a;) = 9{x) 

Jo 2dx J 

or cquivalently 

| [ e i6 ^dy\ 2 = x 2 , x > 0. 



Jo 

But this is only possible if 9{y) is constant, w.l.o.g. we may assume 9 = i.e. u(x) is nonnegative, 
which in turn yields p = since u ^ by assumption. 

3.3. Summary. Let us summarize what we have shown in this section. 

Theorem 3.1. Assume the disease case io = p 2 /(3, a = po/(3. The the semigroup T(t) = e~ Lt 
is strongly ergodic, it converges strongly to the projection V onto the kernel N(L) of L along 
its range R(L). The kernel is one- dimensional and spanned by e(x) = {(3 / po) 2 $>(/3x / po) , where 
$(z) = (z + z 2 /2)e~( z+z ' 2 ' , and the projection V is given by 

/>oo 

Tu(x) = [ / (a + y)u(y)dy]e(x) =< e*\u > e(x), x > 0, u e X. 
Jo 

Remark. We do not know whether the ergodicity is exponential since it is not clear that the 
type of the semigroup e~ Lt restricted to R(L) is negative. 

4. Well-posedness of the Non-Autonomous Evolution 

4.1. The Trivial Evolution. Let lj € C(R+) be positive, such that < uj^ = lim.t_ ) . 00 uj{t) 
exists, and assume uj(-) — Woo G ii(M+). Let 

u>+ = maxu(s) and w_ = minw(s), 

s>0 s>0 

and note that u + > uj_ > 0. We are particularly interested in the cases = Ar/7, the disease- 
free case, and lo^ = Pq/ (3, the disease case. We want to show that the nonautonomous problem is 
well-posed in X = ii(R + ; (a + x)dx). We begin with the problem 

(4.1) d t u(t, x) + uj(t)d x u(t,x) + (p + (3x)u(t, x) = 0, x > 0, t > s > 

u(s, x) = g(x), u(t, 0) = 0, t > s > 0, x > 0. 

The method of characteristics yields easily the evolution operator Uo(t,s) for this problem. It is 
given by 

(4.2) [U (t, s)g] (x) = u(t, x) = g(x - f Lj(r)dr)e-^ s ^ , 

J S 

(j)(t,s,x) = p (t- s)+/3{t -s){x- J Lo(T)dr) + (3 J (t - r)w(r)dT, 

if we extend g trivially to R. We obviously have the estimate \Uo(t, s)|g(x) < e -t*o{t-s) ^ an j u ^ x ^ 
is a strong solution in X if the initial function g belongs to D defined by 

D := {g e Li(R+) : x 2 g,g',xg' S i 1 (R+),. 9 (0) = 0}. 

We also need the solution of 

d t u(t, x) + ui(t)d x u(t, x) + (/i + f3x)u(t, x) = 0, x > 0,t > s >0 

(4.3) u(s,x)=0, u(t,0) = h(t), t> s>0,x>0. 
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Again the method of characteristics applies and yields with K(t, x) — f, t *(r — p(t, x))uj(r)dr the 
formula 

[V {t,s)h](x)=u(t,x) = hipit^^e-^-^^+^-^^-^^^l 
for x < ui(r)dr, and zero elsewhere, where the function p(t,x) is defined by the equation 

(4.4) x = f uj(r)dr; 

Jp 

note that this equation has a unique solution p(t, x) 6 (s, t), since uj(t) > cj_ > for all r > 0, by 
assumption, and x < J to(r)dr. Observe that with Ko(t,s) — J s u(r)dr we have 

3 l-Ko(t.s) 

(a + x)[V (t,s)h](x)dx < / (a + x)e-^ t -P ( - t ' x ^dx 



<\h\oo (a+ uj^dryr^-^Luipit.x^do 



< |^| 00 w + / (a + u+^e-^do < CI^U, 
Jo 

by the variable transformation a = p(t,x). Thus the part coming from a nontrivial bounded 
boundary value h is bounded in X . 

4.2. Well-posedness for the Full Problem. Let us now consider the full problem, i.e. 

poo 

(4.5) d t u(t, x) + u(t)d x u(t, x) + (p + /3x)u(t, x) =2/3 u(t, y)dy, 

J X 

u(s,x)=g(x), u(t,0) = 0, t>s>0,x>0. 

Since the standard cone in X is reproducing, i.e. L\ = L\ — L^, we may restrict attention to 
nonnegative initial functions g. We define the sequence u n inductively by 

«i(t) := U (t, s)g, u n+1 {t) = Ul (t) + / U (t, r)Bu n (r)dr, t>s>0. 



Since Uo(t,s) is positive the functions u n are as well, and U2(t) > Ui(t) since B is positive. 
Inductively we obtain with 

u n+1 (t) - u n (t) = / U (t,r)B(u n (r) - w„_i(r))rfr, t > s > 0, 

J s 

that the functions u n are pointwise increasing w.r.t. n G N. 
Suppose that g € D. Then u n is a strong solution of 

f 00 

d t u n (t,x) + Lj(t)d x u n (t, x) + (po + (3x)u n (t,x) — 2(3 / u n -i(t,y)dy 



<2(3 u n {t,y)dy, x>0,t> s>0 

J X 

u(s, x) = g(x), u(t, 0) = 0, t > s > 0, x > 0, 

i.e. u„ is a strong lower solution of (|4.5|l . Multiplying the equation with cc 1 and integrating over 
K + this yields with Zi(t) = \x l u n (i)\i 

d t z (t) + p z (t) + 0zi(t) < 2p Zl (t), 

for i = 0, and for i = 1 

d t zx(i) - w(t)2b(t) + ^i (t) + /M*) < Pz 2 {t). 

Setting z(t) = (z (t), zi(t)) T , b(t) = (0, — cjoo)zo(i)) T , and defining G by the 2 x 2-matrix 
with entries — fio, 0,u!cx>, — Mo> this inequality becomes 

d t z(t) < Gz(t) + b(t), t>s>0. 
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The eigenvalues of G are given by A± = — (Xq ± \/ /3ll>oo which are both nonpositive if 0u> ^ < ^q, 
which is true in both, the disease-free and the disease case. Since e Gt is positive we may conclude 

z{t) < e G(t - s) z{s) + f e G ^b{r)dr. 

J s 

Boundedness of e Gt then implies an inequality of the form 

\z(t)\<C + C [ \w(r) - uJocWzir^dr, t>s>0, 



which implies boundedness of z(t) on [s, oo) since (w(-) —lu^) € Li(R+) by assumption. Note that 
the constant C depends only on the parameters (J>o,fl, k>oo and on \ \g\\. 

Therefore the functions u n (t) are bounded in X uniformly in t and n. By monotone convergence 
we may conclude u„(t) — ► u(t) in X for each t > s. Since £? is positive, Bu n — > Bu in Li(K+) as 
well, and then also 

(4.6) tt(t) = C/ (t, s)p + / C/ (t, r)Bu(r)dr, t>s>0, 



at least in Li(R + ). A density argument finally shows that this conclusion is valid for all initial 
data g E X. 

Remark. It is not clear that solutions of 1)4.6(1 are unique. The reason for this is that B is 
unbounded. Therefore we need another definition of mild solution. 

Definition. Let / G Lij oc (R + ; X) . 

(i) We call a function u € C(R+; X) strong solution of 

d t u(t,x) +uj(t)d x u(t,x) + (fjLQ +{}%)u(t i x) = 2f3 u(t,y)dy + f(t,x), 

J X 

(4.7) u(s,x) = g(x), u(t,0) = 0, t> s>0,x > 0. 

i/u£ C 1 (K + ; X) n C(M + ; D) and \4.7jj is valid pointwise. 

(ii) We call a function u £ C(M + ;X) mild solution of |^.7| ) if there are f n £ Lu oc (R + ; X) and 
strong solutions u n of |^.7| ) such that u n — > u anrf f n —*f asfi->oo, m X, uniformly on compact 
intervals. 

Suppose that g E D has compact support. Then each iteration u n (t) has also compact support, 
namely 

supp u n (t) C supp g + u + [0, <] , 
for each n E N. Therefore each function u n (t) is a strong solution of l|4.7|l with inhomogeneity 
f n {t) — B{u n -\{t) — u n {t)) ■ This proves that the limit u{t) is a mild solution. Approximation then 
shows that (|4.5() has at least one mild solution, for each initial value g S X . 

Uniqueness of mild solutions can be obtained as follows. If u is a strong solution of 14.7JI then 
the equation yields as above the inequality 

ft||«(t)||< W+ ||«(t)|| + ||/(i)||, t>0, 

hence 

|K*)||<e"+(*- s )||< ? ||+ f e^-^WfirWr. 



By approximation this inequality is also valid for mild solutions, hence u = in case / = g = 0. 
Thus mild solutions are unique and of course they satisfy the integral equation l|4.6|l . 

4.3. Summary. We have proved the following result about well-posedness of (14.5(1 

Theorem 4.1. Suppose to 6 C(R+) is a given strictly positive function, such that lOqo = 
limt^oo uj(t) > exists and w(-) — lOqo € £i(R+). Then (14.5(1 is well-posed in the sense of the 
definition given above. There exists a unique evolution operator U(t,s) in X generated by 1(4.5(1 . 
which is bounded in X , uniformly in < s < t < oo, and positive. Moreover, (14.5(1 has finite speed 
of propagation with maximum speed less than oj + = sup t>0 Lo(t). 
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4.4. Higher Order Bounds. Consider an initial function g £ C^°(0,oo). Then u\ is smooth as 
well and has compact support for each t > s. Then the same holds true for U2, hence by induction 
for all u n . Setting v n = d x u n we have the following problem for v n . 

(4.8) d t v n + uj(t)d x v n + (/j, + (3x)v n = ~P[u n + 2u„_i], 

v n (s,x) = g'(x), v n (t, 0) = ip n (t), t>s>0,x>0 

where ip n (t) = ^|u„_i(i)|i. This implies 

d x u n {t) = v n (t) = U (t, s)g' -pj U (t, r)[u n (r) + 2u„_ 1 (r)]dr + w n (t), t>s>0, 

with 

w n (t) = 2f3V (t,s)[\u n - 1 (-)\i/u(-)}. 

Uniform boundedness of u n in X and exponential stability of the evolution operator Uo(t, s) in X 
then implies boundedness of d x u n in X. Passing to the limit we get 

d x u(t) = U Q (t, s)g' - 3/3 [ U (t,r)u{r)dr + w(t), t>s>0, 

J S 

where 

w(t,x) = 20V o (t,s)[\u(-)\ 1 /uj(-)]. 

This yields d x u £ Cb([s, oo); X). The last identity was proven for g £ C^°(0, oo), but via density 
can be extended to g £ D. 

5. Convergence 

We are now ready to prove the main result on convergence. Let us first look at the disease- 
free case. Then with A(t), B, defined as in section 2, and Lit) — A(t) — B, we know that L(t) 
is strictly accretive for large times t if the parameter a is chosen in a £ (Xt/jhq, iaq/0). This 
proves exponential stability of the trivial solution in the disease-free case, with decay rate at least 
Mo - y/X(3r/j. 

Suppose we have a solution u of the nonautonomous problem in the disease case such that d x u(t) 
is bounded in X . Then we may write 



(5.1) d t u + Uood x u + (fi + Px)u - 2/3 / u(t,y)dy = (wqc - w(t))d x u, 

J X 

u(0, x) = g(x), u(t, 0) = 0, t > 0, x > 0. 
Therefore we obtain the identity 

u(t)=e- Lt g+ [ e- £ (*- r '(woo - Lu{r))d x u(r))dr, t>Q. 



JO 

We know from Section 9 that e~ Lt converges strongly in X to the ergodic projection V . On the 
other hand, the scalar function lu(-) — w M belongs to Li(R_|_) by assumption. This then implies 

u(t) ->Moo G R{V). 

Thus we have convergence in X to a unique element for all nonnegative solutions with initial values 
in D. Since the evolution operator associated with (|4.5() is bounded in X , this convergence extends 
to all initial values uq £ X . 

Returning now to the system we may compute the limit Moo. For this purpose recall that 

U(t) = J xo u{t,x)dx — > Uoo and P(t) = J xo u(t,x)xdx — > P^. This implies 

Moo = Urn Vu{t) = lim [aU(t) + P(t) - x U(t)]e = [fiU^/P + P^je. 

t — >oc t — >oo 

Note that Uoo is independent of the initial values Vq and uq. 
This completes the proof of Theorem 1.2. 
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